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Abstract 


Recent ideas about understanding physics of black hole information¬ 
processing in terms of quantum criticality allow us to implement black 
hole mechanisms of quantum computing within critical Bose-Einstein 
systems. The generic feature, uncovered both by analytic and nu¬ 
meric studies, is the emergence at the critical point of gapless weakly- 
interacting modes, which act as qubits for information-storage at a 
very low energy cost. These modes can be effectively described in 
terms of either Bogoliubov or Goldstone degrees of freedom. The 
ground-state at the critical point is maximally entangled and far from 
being classical. We confirm this near-critical behavior by a new an¬ 
alytic method. We compute growth of entanglement and show its 
consistency with black hole type behavior. On the other hand, in the 
over-critical regime the system develops a Lyapunov exponent and 
scrambles quantum information very fast. By, manipulating the sys¬ 
tem parameters externally, we can put it in and out of various regimes 
and in this way control the sequence of information storage and pro¬ 
cessing. By using gapless Bogoliubov modes as control qubits we de¬ 
sign some simple logic gates. 
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1 Introduction 


The recent advances [l]^5] in understanding physics of black hole information 
processing in terms of the phenomenon of quantum criticality opens up a 
possibility of studying and implementing its key mechanisms in ordinary 
Bose-Einstein systems, such as for example, cold atomic gases. 

Black holes are fascinating objects from the point of view of quantum infor¬ 
mation processing. The information storage capacity of a black hole of radius 
R is measured by Bekenstein entropy |6j, which scales as area in Planck units 


N 


R 2 

U p - 


( 1 ) 


The Planck area in terms of Planck and Newton constants is defined as 
Lp = HG N . Thus, the black hole entropy diverges in the classical limit, 
h —> 0, if R and Gn are kept finite. In this limit black hole half life-time, 
tbh ~ iV 3//2 Lp, also becomes infinite. Hence, a classical black hole can store 
an infinite amount of information eternally. 

What sort of microscopic physics is responsible for this type of behavior? 
According to [I], the internal qubit degrees of freedom that store information 
in a black hole satisfy the following properties. The energy gap in each qubit, 
is suppressed as some positive power of 1/N relative to “naive” energy gap 
that would be expected for an ordinary weakly-interacting quantum system 
of a similar size (for example, for a gas of free cold bosons of mass m in a 
box of size R , the gap between the levels would scale as ~ h 2 /2mR 2 ). In 
the same time the interaction strength of the qubit degree of freedom, with 
itself and the rest of the system, also must be suppressed by powers of 1/N, 
in order to accommodate the fact that in N —► oo limit the state of the 
qubit is an eigenstate of the Hamiltonian and the information stored in it is 
frozen eternally. It was also suggested, that such a behavior emerges from the 
following microscopic physics. The black hole can be described as a bound- 
state of N soft gravitons at the quantum critical point |1|. The quantum 
critical point amounts to the fact that the graviton-graviton coupling a and 
their occupation number N satisfy aN = 1. That is, the critical point 
describes a quantum transition between the free graviton Bose-gas and a 
collective bound-state. It is this quantum criticality that gives rise to almost- 
gapless collective excitations that play the role of the information qubits. [^] 

2 Throughout the paper we use the term “qubit” in generalized sense, although the near¬ 
gapless excitations living at the quantum critical point are in general multi-level quantum 
systems. This allows for more possibilities of quantum information storage in these modes 
(e.g., using their coherent states 0). but, for simplicity, we can restrict ourselves to two 
basis-states per mode that can be effectively separated from the rest. 
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We shall now abandon the gravitational aspect of this picture and take its 
main message, that the quantum criticality of Bose-gas is the key to the un¬ 
usual features of black hole information processing. From here it follows that 
similar properties must be exhibited by much simpler systems that also can 
be placed at the quantum critical point, such as the critical Bose-Einstein 
condensates of ordinary particles. This expectation was confirmed by the pre¬ 
vious studies, [1-3,5], where it was shown that already the simplest prototype 
models, indeed exhibit some of the key features of black hole informatics. 

One of the generic features is the appearance of nearly-gapless weakly-interacting 
qubits. The important property is that the energy gap of qubits is controlled 
by closeness to criticality and can be made arbitrarily small. In the same 
time, interaction of qubits is suppressed by 1/N and can be made arbitrarily 
weak by choosing system with large N. 

Another generic feature is the high level of entanglement at the quantum 
critical point |2|. It was also shown |3] that in the over-critical phase, where 
a Lyapunov exponent and chaos develops, the system becomes a very fast 
scrambler of quantum information This result gives an explicit microscopic 
justification to the conjecture of 10 about the black hole’s fast scrambling 
property. It also shows that this property is not limited to black holes, but 
is rather generic for critical Bose-Einstein systems with Lyapunov exponent 
and a chaotic behavior. 


All the above studies show that the effective theory of information storage 
near the quantum critical point is represented by gapless and very weakly 
interacting qubit degrees of freedom. In these qubits one can store quantum 
information cheaply and for a macroscopically-long time. 

In the present paper we take further steps in establishing the role of quantum 
criticality in black hole type information processing. In particular, we focus 
on implementing this way of quantum computing in simple Bose-Einstein 
systems. 

We will be interested in the regime in which a system is very close to the quan¬ 
tum critical point. That is, in a situation when a properly-normalized collec¬ 
tive coupling, A = aN, is extremely close to one. As said above, this regime 
is universally characterized by appearance of almost gapless qubits. For ex¬ 
ample, in Bogoluibov approximation the energy gap scales as e = \/l — aN. 
However, the appearance of the nearly-gapless mode is well-established also 
beyond this approximation, both by numeric |2,3,5] as well as analytic meth¬ 
ods [5]. When approaching the critical point from the weak coupling (e —*■ 0+) 
the system is characterized by a slow evolution of qubits and a very sharp 
increase of the entanglement of the ground-state at A = 1. Below we shall 
re-establish this result analytically by using a new method. 
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In the overcritical regime, e becomes imaginary. There the system develops 
a Lyapunov exponent and becomes a fast scrambler. 

By manipulating the parameters externally we can adjust proximity to the 
critical point and take the system into different regimes. In this way we can 
control the efficiency of information-storage and processing. 

We construct an explicit sequence of information storage and read-out in 
near-gapless Bogoliubov modes of Bose-Einstein condensate at the quantum 
criticality. Using these modes as control qubits, we construct some quantum 
logic gates. We outline an explicit evolution sequence that enables to store 
information in these qubits and later decode it. 

We would like to note that the prototype Bose-Einstein systems, can be 
highly useful for exploring other aspects of the black hole’s quantum-critical 
portrait, such as, e.g., collapse and evaporation, 0 . § In the present pa¬ 
per, however, we shall focus solely on the physics of quantum information¬ 
processing. 


2 Cheap Qubits at Near Criticality 

We start with the system that was used in jT] as the simplest prototype 
model for exploring the connection between black hole information processing 
and quantum criticality. It represents a gas of cold bosons with a delta- 
function-type attractive interaction. The Hamiltonian in d-dimensions has 
the following form, 


n 



^h 2 a 

2m 


i> - gh 


J d d Xlp + lp + , 


( 2 ) 


where ip = ^] k ~^e z R x a*,, V = ( 2nR) d is the d-dimensional volume and k is 

the d-dimensional wave-number vector. are creation and annihilation 

operators of bosons of momentum-number vector k. These operators satisfy 
the usual algebra: [a*,, a) k P\ = dkk' an d all other commutators zero. The 
parameter g controls the strength of the coupling. The above system exhibits 
a quantum critical behavior. In particular, in d = 1, it is known to exhibit a 
quantum phase transition towards the bright soliton phase, studied in details 
in |9l. 

It is useful to represent the Hamiltonian in the form, ~H = 9/ ^ m H. The 
quantity 2 R 2 m i s a convenient unit for energy-measurement, since it gives us 
a clear idea about the energy-cost of a given process, relative to the kinetic 
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energy of a first non-zero momentum mode of a single free boson. For this 
reason all our further discussions will take place in the units 2 j&rn = 1- That 
is, we effectively switch to the Hamiltonian H. 

Introducing a notation a = (y^) and taking d = 1 this Hamiltonian 
takes the form, 

H = Yjk 2 a\(i k -j 2 a\a\a k3 a k i. (3) 

k ki+k2~k3—k4=0 


We restrict ourselves to three levels, k = 0, +1. Legitimacy of this approxima¬ 
tion for the regime of interest was firmly established from the previous analy¬ 
sis of this system Bill TT]. Let us define a triplet operator a; = (a_i, a 0 , ai). 
Correspondingly, we define the number operators n t = ajaj. Since the parti¬ 
cle number is conserved we can restrict ourselves to an iV-particle sector. 
Then the Hamiltonian in terms of the triplet a; can be written as [5], 


H = 




n,; 


^ (aja^cioao + a\a} 0 aia_i) + j ^ nf + H(ri ). (4) 


i =—1 


where, H{n ) is the part of the Hamiltonian that depends on the total number 
operator n = n_i + uq + n\ and has the following form, 


H(n) 


a n a . 

- n + —n + n in 

2 4 


N), 


(5) 


where /x is a Lagrange multiplier. For aN « 1 the minimum of the Hamilto¬ 
nian with a very good approximation is achieved at the state n, = (0, N, 0). 
Expanding around this state the unsuppressed part of the Hamiltonian can 
be expressed as the mass matrix that mixes the operator doublets (a},a_i) 
and (ai, al-J, 

( 6 ) 


(l - %N, -f N 
{-‘IN, 1-f N 


We shall now diagonalize this bilinear hamiltonian via a Bogoliubov canonical 
transformation, 


where, 


with 


'+1 = ua ±l - va '+l , 

(7) 

1 + e 

1 — e 

(8) 

2Ve ’ 

CM 

II 

e = Vl — cxN. 

(9) 
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We thus arrive at the following effective low energy Hamiltonian for Bogoli- 
ubov modes: 


H b - €((-1,6-, + b\M + 0(b'), 


( 10 ) 


where 0(6 4 ) stands for momentum-conserving quartic interactions among 
6 +1 and b\ x operators. 

The proximity of the system to the quantum critical point can be measured 
by the parameter e. This parameter can be tuned to be arbitrarily small by 
choosing aN sufficiently close to 1. Moreover, for any given value of e > 0 
we can always take large-enough N so that the interaction term in (10) can 
be neglected for low-enough 6-levels. We shall refer to this choice of system 
parameters e and N, as a double-scaling regime of the Bogoliubov domain. 
The validity range of this regime will be quantified below. Before doing 
this, we would like to pause and reflect on the fundamental importance of 
the Hamiltonian Hb for quantum information storage and processing. This 
connection was already much appreciated in [l]-3j|5] and is also central to 
the present work. Namely, in the double-scaling regime Hi, becomes a free 
Hamiltonian of 6-modes with energy gap e. Thus, the states of small-enough 
occupation number of 6-particles, \n)b, are approximate energy eigenstates, 
Hb\n)b = E n \ri)b of eigenvalues E n = ne. Here and below n denotes an 
eigenvalue of the 6-number operator, 6 + 6. As said above, for each value of 
e this approximation can be made arbitrarily accurate by taking N to be 
sufficiently large. The precise scaling will be quantified below. 


Thus, using 6-modes as qubits, the energy-cost of information-storage can 
be made arbitrarily cheap as compared to typical energy gap in the ordinary 
non-critical systems. Secondly, the information stored in 6-modes stays intact 
for a macroscopically-long time, controlled by the powers of N. 


3 Regimes 

Let us now quantify the parameter-choices for various regimes. We start 
with the double-scaling regime of Bogoliubov domain in which, as described 
above, Bogoliubov effective theory can be made applicable by taking both e 
and 1/N sufficiently small. We shall identify the required relation between e 
and N. 

In order to achieve this, let us notice that for any fixed values of N and e 
there exists a maximal occupation number n max , above which the Bogoliubov 
treatment breaks down. Obviously this value of n also measures the number 
of states that fall within Bogoliubov treatment. This number can be easily 
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estimated from the following argument. 

First, we shall require that the energy of the state must be within the gap 
A Hi, ~ 1. From the form (10) it is clear that this condition translates as the 


bound on occupation number of 6-modes, 


n < e 


-i 


(ii) 


Secondly, the validity of Bogoluibov approximation imposes a condition that 
the interaction term among 6-modes must be sub-dominant with respect to 
the first term in (10). This condition translates as a second bound on n, 


n < e 3 N. 


( 12 ) 


The maximal value n m(7T that satisfies both conditions is bounded from above 


as 


n r , 


< 


N ~ l/4 . 


(13) 


Thus, for any N we can adjust e ~ 1V~ 1//4 in such a way that there are 
~ A 1 / 4 approximate energy eigenstates within the A Hj, ~ 1 energy gap. The 
storage of information within the lowest lying states is extremely cheap and 
the evolution time is long. 

We can determine the absolute bound on the domain of Bogoliubov treat¬ 
ment, by requiring that there exist at least one number-eigenstate within its 
validity. That is, we take n ~ 1. Then from (12) we get the absolute lower 
bound on e, 

e > A" 1/3 . (14) 

Thus, by taking N sufficiently large we can make the energy gap arbitrarily- 


small and controllable within the effective theory given by (10). 


The bound (14) marks the domain of applicability of the Bogoliubov treat¬ 


ment. However, near-gapless states continue to exist even beyond this do¬ 
main. In fact, the alternative methods [3,|5] as well as the new one that will 
be discussed below show that density of such states in the vicinity of critical 
point remains huge. However, they cannot be analyzed within the Bogoli- 
ubov effective theory given by (10). In this domain e is still real, but it is 
below the Bogoliubov bound, e < N ~k 3 . We shall refer to this domain as to 
trans-Bogoliubov. In this domain the Bogoliubov approximation is no longer 
valid and we have to employ a different method. In |2,3| this was achieved by 
exact numerical diagonalization and in [5] by developing an effective Nambu- 
Goldstone description of modes. Below we shall use a different method for 
treating this domain. 
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Finally, there is a third domain, aN > 1, in which e becomes imaginary. In 
this domain, the state in which only no-level is macroscopically-occupied be¬ 
comes unstable and system relaxes to a new ground-state of a bright soliton. 
As shown in p], during this instability the system is a very efficient scram¬ 
bler of information due to its Lyapunov exponent and a very high density of 
states near the critical point. 


To summarize, we identify three domains of interest in the vicinity of quan¬ 
tum criticality. The two belong to the sub-critical regime (aN < 1, i.e., 
e > 0). The first of these two is the double-scaling Bogoliubov domain, in 
which e is real and satisfies the bound (14). The second is trans-Bogoliubov 
domain, in which 0 < e < N _1 / 3 . Both of these domains are characterized 
by a very low energy gap as well as by slow evolution of low-lying states. In 
particular, generation of entanglement takes macroscopic time, enhanced by 
powers of N. 

Finally, there is an over-critical fast-scrambling domain with imaginary e, in 
which the system becomes unstable. Unlike the other two this domain is 
characterized by an extremely fast growth of the entanglement and scram¬ 
bling of information. 


Having identified the different domains, we can manipulate the system putting 
it in and out of different regimes in accordance with the information-processing 
task. This can be done, because by manipulating characteristics of the Bose- 
Einstein system we can change criticality parameters, such as e (or aN). 

For example, consider the following sequence of information storage and pro¬ 
cessing. First, keeping system within Bogoliubov domain, we can dial 6- 
modes by coupling the system temporarily to some external influence. This 
dial-up encodes the initial message in form of the occupation numbers of 
6-modes. Then, after switching-off the external influence, information re¬ 
mains stored in 6-modes for a long time due to very slow evolution. Next, 
by switching on the coupling with some external qubits, we can process the 
stored information by using 6-modes as control qubits. We give examples of 
designing some simple logic gates in the next section. 

Alternatively, putting the system in over-critical domain, we can scramble 
the information and transform the original message in a highly-entangled 
state. In this regime e becomes imaginary. The system acquires a Lyapunov 
exponent and scrambling time becomes logarithmic in N. In principle, by 
manipulating the parameter e externally, one can control the level of entan¬ 
glement. 



4 Qubit Evolution 


As an example of using Bogoliubov modes as qubits near the critical point, 
consider the following system. Let the Bogoliubov mode operator be 6. Near 
the critical point in the double scaling limit the Hamiltonian (10) can be 
approximated as a free Hamiltonian of 6-modes. Let us consider a single 
lowest lying 6-mode, 

H b = e 6^6 . (15) 


We assume that V -1 / 4 > e > 1/3 in order to be within the validity 

domain of the Bogoliubov approximation and ignore 1/V-suppressed inter¬ 
action terms. 

Let us couple 6 to a probe external quibit c with creation annihilation algebra 
[c, c'] = 1. We choose the Hamiltonian in the following form, 


H = eb^b + b^b{yc + z/*(d) + 2 Sc^c. (16) 


Here 5 is a real parameter and v is complex, which for simplicity we choose 
to be real as well. For non-zero coupling v the two degrees of freedom can 
control the time-evolution of each other. In particular c can be used for 
encoding information in a state of 6-modes, as well as for reading-out this 
information. 


4.1 6-Mode as Control Qubit 

Let us first consider a situation in which 6 acts as a control qubit for c. For 
this purpose it is enough to take c as a two-level system described by the 
states |0) c and |l) c . 

For \v\ « y/e~5, the ground-states of the two systems are uncorrelated and 
are given by a tensor product of two free-oscillator ground-states, |0)b Cg) |0) c . 
However, for \v\ » Ve~6 the two systems become strongly correlated. 

In particular, in the regime e ~ v » S, the ground-state of the qubit c is 
determined by the state of 6. Thus, 6 acts as a control qubit for c. 


It is easy to compute the time evolution as given by (16). If we label the 
basis states as 

a 


A := a | m) b | 0) c + f3 \ m) b \ l) c , 


we find 



(17) 
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(19) 

For time dependent v = v(t) one merely has to substitute tv with $ v{t)dt. 
We see that the Hamiltonian acts as a generalized controlled gate, in the 
sense that there are two different time scales for qubit state with 0 occupa¬ 
tion number of 6-particles (m = 0) and with one or more 6-particles ( m > 0). 
For the m = 0 state the time scale is | whereas for m > 0, after the adiabatic 
turning on of v(t ) until the final value v en d » 5, the time scale is given by 
, 1 ss —-— « f. For v en d » 5 the m = 0 state can be viewed as 

(iTl t'end) 2 +<5 2 m ^ 5 

stationary. 

For example, restricting to two levels m = 0,1 in the limit 6 = 0 and for 
tv = this system acts as a standard CNOT gate, because it acts trivially 
on the state with m = 0, whereas for m = 1 we have (up to a phase factor), 



4.2 Encoding Information in 6-Mode via Time Evolu¬ 
tion with Coherent States 

For encoding the initial information in 6-mode an external mode with similar 
coupling can be used. We keep the notation c for this mode. We shall 
consider the case when the state of 6-mode is controlled by the coherent 
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state of c. Preparing the state of c in a coherent state \^) c and switching on 
coupling v(t), we can evolve the 6 -mode to a particular state controlled by 
the coherent state parameter 7 . We take the Hamiltonian to be (16) again, 
however, now the operator c is no longer restricted to two levels, but is a 
fully-fledged harmonic oscillator. We want to study the time evolution of 
6 -particle states tensored with coherent states of c. Again the problem can 
be treated analytically in simple terms. We now compute 


-i/ft 


m )b ® I l)c . 


( 20 ) 


where | y) c is a normalized coherent state of c, i.e., c | 'y) c = 7 | 7 ) c . We get 
(H = | m)j ® e -“< m ‘'( c+c, )+Mct c ) | ^ 


( 21 ) 


We now perform the following trick. Let us define d := c + Note that 
coherent states for c are also coherent states for d, namely 


d I 7> c = { 7 + ijfj I t) c <*=* I 7> c = 


mv 

7 + ^ 


Then 


and we hnd 


T 777 

25c f c + mv(c + c*) = 25d f (i- - 

2 o 


if(mi/(c+ct)+2<5ctc) | ^ ^—it28dftd 


mv 
7+ 25 


( 22 ) 


(23) 


(24) 


The last equation is simply the time evolution of a coherent state with fre¬ 
quency 25, so that in total we hnd: 


-i Ht 


m )b®\l)c = e 


_ -i t(em- 


(mu) 

25 


m) b ® 




(25) 


The time evolution preserves coherence and the speed of the evolution de¬ 
pends on [7 + ^|. Incidentally we observe that states 7 = — 77 do not 

evolve, i.e., are energy eigenstates with energy (em — ^s~)- H 


We note that we can dial 6 -states, e.g., by first preparing a particular coher¬ 
ent state in c, say Q]/7j | m j) b ) ® | — z ^) , where /3j are some parameters 

3 Obviously, the Hamiltonian is unbounded from below for time-independent v 
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describing the state of b, and then after some time probing it against itself. 
The matrix element 


m 0 u 


® b {X 



m j) b ® 



(26) 


is dominated by the contribution from | mo) b for any | X) b (as long as it has 
a non-zero overlap with | mo)J since 


m qv 


-U2S (rrij-m^v 

! 2h 


m 3 U \ = e -^K-^o) 2 sin 2 (5i)-i(...) 
25 /c 


(27) 


is exponentially suppressed for v » 5 for any rrij ^ itiq when averaged over 
time. Here, the phase factor in the exponent is (...) = mo (mo—mj)sm(2t8). 
Thus, if one has a pointer which reacts to classical configurations of c, 
by preparing and measuring states with such a pointer one can collapse 
randomly-generated states to fixed b mode states - i.e., one can prepare cer¬ 
tain 6 -mode states by performing measurements on c only. 


5 Beyond the Bogoliubov Domain 

We now wish to study evolution of the system near the critical point in 
the domain 0 < e < iV -1 / 3 , which is beyond the validity of the Bogoliubov 
double-scaling domain. In |2,3,B] this regime was studied numerically. In (5) 
a new method was developed in which the quantum phase transition was 
mapped onto a Goldstone phenomenon in a sigma-model that describes a 
spontaneous breaking of the global U (3)-symmetry of Hamiltonian (|B]), under 
which the triplet of operators (a_i,ao,cq) transforms in the fundamental 
representation. In the state (0, N, 0) this symmetry is spontaneously broken 
down to 1/(2). The quantum phase transition corresponds to a change of 
symmetry breaking pattern to a smaller symmetry group achieved on the 
ground-state (x, N—2x, x). As shown in [B] the transition to this new ground- 
state with x 7 ^ 0 takes place beyond aN = 1 point. This mapping allows 
to relate the gapless mode, appearing at the quantum critical point, to a 
Nambu-Goldstone mode of the sigma-model that becomes massless at the 

same transition point. The effective Hamiltonian of the gapless mode up 

n 2 

to 1/iV-corrections is H go id = where n go id is the occupation number 

of the gapless mode. The consistency of this effective description was also 
checked by numerical study of various characteristics, such as, derivation of 
the ground-state, calculation of the energy for different occupation numbers 
of Goldstone mode and exact time-evolution of states [B]. 


12 











Thus, the effective Goldstone approach shows that the existence of nearly- 
gapless weakly-coupled excitations near the critical point persists both in 
trans-Bogoliubov domain, 0 < e < N^ 1 ' 3 , as well as in the over-critical one, 
as long as |e| (equivalently x) is small. 

Below, we shall complete the study by applying a new method, which allows 
to treat the system analytically beyond the Bogoliubov approximation in 
sub-critical domain and compute certain important characteristics, such as 
the growth of entropy at the critical point, explicitly. In this paper we shall 
only outline the technique that is essential for our analysis. The detailed 


description of the original method will be given in 11 


5.1 The Setup 

As discussed above, after we restrict ourselves to an IV-particlc sector with 
total momentum zero and cut off the modes with momenta \k\ > 1, we arrive 
at the Hamiltonian Q. The Hilbert space becomes finite dimensional. Its 
basis can be written as 


| ?z) : = | n_i = n, no = N — 2n, ni = n), 

where n goes from zero to y. The three mode Hamiltonian in this basis is a 
(y + 1) x (y + 1) dimensional tridiagonal matrix of the following form: 


where the entries are 


H = 


( do 
hi 
0 


hi 

di 

h 2 


0 ...\ 



-cN, 



r / ,, 3A n 2 

d n = (2-A )n+ y — 


, A A 2 

hn = ~2 n+ N n ' 


(28) 


(29) 


Here A = aN is the same parameter as defined in the introduction and c is a 
positive constant Q 

The spectrum and the eigenstates of H were computed analytically by means of 
a novel diagonalization technique. We will not explain this technique here, but 
rather state the dictionary provided in |TT] shortly and refer the reader to the 
original paper for more details. This dictionary is 


n) <—» b n (x ) := e 2 L n (x) 

H — - 1) + 2(A - 1) ( xdl + d x + 2 “ X ') + ±H int . 


4 We have approximated y/(N — 2 n){N — 2n — 1) « (N — 2 n). 
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5.2 The critical entropy growth 

5.2.1 Entanglement of the condensate state 

At the critical coupling A = 1 the unsuppressed part of H corresponds to a multi¬ 
plication operator. The suppressed interaction part jjH int is extremely important 
for the computation of the eigenvalues and eigenstates of H - indeed it controls 
the depletion and regulates the phase transition. However, for states in which only 
low | n) are occupied the interaction part does not contribute to the time evolution 
for long times for large Therefore we can use the multiplication operator for 
computing the time evolution of those states. This greatly simplifies the problem. 

We wish to compute the time evolution of the condensate state | 0), in other words, 
the coefficient 

K n '■= (n\ e~ lHt | 0). (30) 

From K n we can calculate the 1-particle density matrix and then the von-Neumann 
entropy of the time-evolved condensate can be calculated as function of time. The 
density matrix p in this context is a 3 x 3 diagonal matrix where the entries are 12 

JV 

Poo = Jj Yj\ku\ 2 {N — 2n) 

72=0 

N 

1 2 

1 V* I 12 

pn = p-i-i = V7 2j l Kn l n - 

iV n =0 


(31) 


Now we use the fact that for a long time-scale (proportional to some positive 
power of N, to be determined below) and for small enough n, we can use the 
multiplication operator to compute K n . It turns out that “small enough” n are 
the dominant contribution to the density matrix, so we are in a very comfortable 
position. From the dictionary above we hence obtain: 

*> * = e ’‘ (2 2 +~tt)L ' ( 32 > 

where the approximation is very good for small n over long time scales. From here 

5 This is justified both by our previous discussion as well as by the existing studies 1,5] 
that show that evolution takes macroscopic time in N. 
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we compute: 


Poo = 


2 N 


N_ 

2 


t 2 


- 1 


P11 = 


4JV 


(33) 


1 - 


V 


t 2 


2IV + t 2 + 4 
f 2 + 4 


The von Neumann entropy is given by — Tr(plogp), i.e., by: 

S = - (poolog(poo) + 2pn log(pn)) . 
For small times and large N we find 


(34) 


Poo 


t z 

2N 


1 


P n 


4 AT’ 


and so for early times the entropy grows as 


S'-log 


VnJ n 


(35) 


Thus the characteristic time for the entropy growth is macroscopic and scales as 
t en t ~ yN. This confirms the earlier hndings 2,3,5 indicating that the evolution 
of lowest lying states remains slow, even beyond the Bogoliubov regime, as long 
as the system is subcritical. It is also in agreement with the derivation of entan¬ 
glement near the critical point given in |2|. 

We compare the analytic expression for S with the one obtained numerically using 


the matrix (28) for N = 250, see the image below. We see that even for such a 


small N the analytical expression for S is valid for long times, i.e., the entropy 
indeed grows until its maximal value with the relevant time-scale being y/N. 

The situation changes dramatically in the over-critical domain, in which e becomes 
imaginary and the system becomes a fast scrambler |3j. 

5.2.2 The modes contributing to entanglement 


It is clear that formula (30) gives an overestimate of the real coefficients, since we 


ignored the depletion regulating interaction terms. Nevertheless, quite remarkably, 
even within this approximation n n does not grow uncontrollably in time. Indeed, 
the maximal value that our approximated K n can attain is 


1 


max(|/c n (i)| ) = 
t n + 1 V n + 1 


n 


,-i 


n + 1 


for large n. 


(36) 


However, since it is n\K n \ what enters into the density matrix we find that all 
modes contribute almost democratically to the entanglement entropy within our 
approximations. 
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Comparison to numerics for N=250 



- Exact numerical result 

— Analytical expression 


Figure 1: Comparison of the analytic expression for S with the one obtained 
numerically using the matrix (28) for N = 250. 


5.2.3 Criticality of the ground state 

In the overcritical regime e becomes imaginary and the homogeneous condensate 
is unstable. The new ground state is a non-homogeneous condensate, described 
as a shift-invariant superposition of bright solitons located at different points on 
a circle [9] |^J The ground state for A > 1 is classical in the sense that for large -N 
it can be well-described by a coherent or a product state for long times. In this 
regime, and close to criticality (A — 1 « 1), the initially-homogeneous condensate 
state | 0) entangles very efficiently, because of the instability, which provides a 
Lyapunov exponent, and because of the high-density of states near the critical 
point [3]. 

We must distinguish this situation from a very high entanglement that takes place 
near the critical point when approaching it from the sub-critical regime (A < 1). 
Although the system can reach a very highly entangled state, the growth of the 
entanglement in this regime is slower than in the scrambling regime and takes over 
a macroscopic time scale ~ y/N. This entanglement does not occur because the 
system’s ground state is another classical state, different from |0). Rather, the 
physical reason for high entanglement is that in the critical regime the ground 
state is maximally far from being classical. In a separate work we will show that 
in the N —> cc limit the n point functions do not factorize for any n and hence 
the ground state cannot be approximated by “classical” (e.g., product, coherent, 
squeezed etc.) states. 

To conclude this section, the glimpses of black hole type information physics are 
clearly visible in the above analysis of the entropy growth near the quantum- 

6 This has interesting implications, like the possibility of preparing Schrodinger’s cat 
states using the gas of bosons. In this paper we shall ignore this subtlety since it does not 
affect the present discussion 
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critical point. From general arguments by Page 1121 it has been expected that 
black hole’s half life time plays the transitional role in information-processing, but 
the microscopic reason behind this has not been identified clearly in the past. The 
meaning that quantum-critical portrait of black hole gives to Page’s time is the 
time of maximal entanglement of black hole constituents, which can be modeled in 
group-theoretic terms as number of steps required to entangle states in the Fock 
space of gapless qubits of the critical point (4]. From the above computation of 
the entanglement entropy evolution, as well as from the previous analysis [2,3,5|, 
it is evident that already the simplest critical condensates capture the properties 
of development of maxima entanglement over the power-law time in N. 


6 Conclusions and Outlook 

From the recent studies 010 the following is becoming more and more evident. 
First, the peculiar properties of black hole information-processing originate from 
the fact that black holes represent multi-particle systems at quantum critical point. 
Secondly, the key properties of black hole information processing are shared by 
other critical systems, such as Bose-Einstein condensates. This fact opens up an 
exciting possibility of implementing black hole type quantum computing in such 
systems, both theoretically and experimentally. 

As it was uncovered by the above studies of the simple prototype models, the 
generic property is the appearance near the quantum critical point of nearly-gapless 
weakly-interacting modes that act as qubits for information storage and processing. 
In different approaches these modes have been effectively described as Bogoliubov 
modes of the cold gas of bosons (e.g., gravitons) |T| or as Nambu-Goldstone modes 
of an isomorphic sigma-model |5|. It is evident that seemingly-mysterious prop¬ 
erties, such as, e.g., cheap information storage for macroscopically-long time and 
fast scrambling of information, are determined by the behavior of these quibits in 
various near-critical regimes. In particular, in under-critical domain, the qubits 
store information very cheaply and for a very long time. The ground-state at the 
critical point is highly entangled, but if we approach the critical point from A < 1 
regime, the development of the maximal entanglement takes a very long time. In 
contrast, in the overcritical regime A > 1, Lyapunov exponent ensures the fast 
scrambling of information |3|. 

All these features, enable us to implement various regimes of black hole quantum 
computing by manipulating the system parameters externally. In this paper we 
achieved the following goals. 

First, by coupling the system to an external influence, we designed a simple com¬ 
putation sequence by first encoding information into the Bogoliubov modes and 
later decoding this information. We showed how the information is encoded due 
to interaction of the system with external c-modes. We gave example of how to 
dial 6-modes to the desired states by performing measurements solely on c-states. 
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This process of inscription can be viewed as an effective description of information- 
encodement in a wide range of critical systems reduced to their bare essentials. 
For example, for black holes, the role of the c-modes can be played by an external 
scattered radiation, whereas for cold atomic systems, this can be a laser light or 
another interacting atomic gas. 

We showed, that by tuning parameters of the system, such as N and e, we can make 
the time evolution of 6-modes arbitrarily slow, and thus, the information-storage 
time in these modes - arbitrarily long. 

Next, using 6-modes as control qubits for c-modes, we can read out the stored 
information. We gave some examples of control logic gates. 

Finally, by using a new analytic method, to be discussed in more details in |11] , we 
went beyond the Bogoliubov approximation and studied growth of entanglement. 
We re-confirm previous results |2,|3j|5], that in sub-critical domain the growth 
remains slow even beyond the Bogoliubov regime. In the same time, the final 
level of entanglement is maximal at the critical point. It is tempting to suggest, 
following |3,4j|, that this maximal entanglement of the ground-state near the crit¬ 
ical point can shed light at the black hole state after its half lifetime (i.e., Page’s 
time [12]). 


The next obvious task would be to realize, the type of quantum computing se¬ 
quence considered here, in an experimental setup. In the view of the progress in 
quantum gas experiments 13 , one could envisage designing the systems that can 
be manipulated near the quantum critical point. Perhaps, this may be attempted 
by suitably adjusting the setups which allow to design coupled cold gases in one di¬ 
mension, such as [l41. Since the type of quantum criticality we considered is based 
on attractive interaction, it would be interesting to see, if the similar quantum 
processing can be implemented in systems with stable magnetic droplets observed 
in [l5], which are known to be associated with the attractive interaction between 
the elementary spin-excitations 116 and nucleate at threshold current |17 . 
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